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Abstract — We describe several techniques for using bulk mat- 
ter for special purpose computation. In each case it is necessary to 
use an evolutionary algorithm to program the substrate on which 
the computation is to take place. In addition, the computation 
comes about as a result of nearest neighbour interactions at the 
nano- micro- and meso-scale. In our first example we describe 
evolving a saw-tooth oscillator in a CMOS substrate. In the sec- 
ond example we demonstrate the evolution of a tone discriminator 
by exploiting the physics of liquid crystals. In the third example 
we outline using a simulated magnetic quantum dot array and 
an evolutionary algorithm to develop a pattern matching circuit. 
Another example we describe exploits the micro-scale physics 
of charge density waves in crystal lattices. We show that vastly 
different resistance values can be achieved and controlled in local 
regions to essentially construct a programmable array of coupled 
micro-scale quasiperiodic oscillators. Lastly we show an example 
where evolutionary algorithms could be used to control density 
modulations, and therefore refractive index modulations, in a 
fluid for optical computing. 



correlations among the spins. This system is computationally 
tractable only if we make certain simplifying assumptions. But 
even then, we cannot compute the exact spatial correlations, 
only the general picture. This example is a particularly in- 
teresting problem because we can compute correlations either 
using detailed quantum mechanics and differential equations, 
or we can utilize automata theory and obtain essentially the 
same result (more on this later). Of course the automata theory 
approach is computationally much faster then the differential 
equation approach, and the real-world process is even faster 
Feynman has implied that the automata theory approach is 
a potentially more realistic description of the dynamics at 
the meso-, micro- and nano-scale, than systems of differential 
equations. The tiny "computational agents" at those scales do 
not compute differential equations. They simply interact with 
their nearest neighbors and swap information, as Feynman says 
[3]: 
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I. Introduction 

There are many physical processes that can be described as 
a computation. For example, crystal growth from nucleation, 
corrosion-dendrites on an electrochemical electrode, a drop of 
ink dispersing in a glass of water, are all physical/chemical 
processes of increasing complexity that can be thought of 
as a computation. Further, since biological systems are part 
of the physical universe, the development of an organism 
from a fertilized egg is also a computational process. The 
common element in each of these processes is the fact that the 
computation is taking place only between nearest neighbors. 
There is no global clock, or central processor to distribute tasks 
to the individual processes comprising the overall system. 

Many of these processes are either difficult or way beyond 
current computational abilities for modelling. Given a super- 
computer, and the set of differential equations and boundary 
conditions that describe some of these processes, it is likely 
that we would find that our computed results are only an 
approximation of the real-world system. The world is a better 
model of itself than the models we can induce from our 
data. Many of these problems are not only computationally 
intractable, but also computationally undecidable [1], [2]. 

As an example, consider an array of magnetic spins in which 
each site takes on only one of two spin states. At a high 
temperature the spins will be randomized, but as we cool the 
array down to much lower temperatures we will find spatial 



"It always bothers me that according to the laws 
[of physics] as we understand them today, it takes 
a computing machine an infinite number of logical 
operations to figure out what goes on in no matter 
how tiny a region of space and no matter how tiny a 
region of time. How can all that be going on in that 
tiny space? Why should it take an infinite amount of 
logic to figure out what one tiny piece of space-time 
is going to do?" 

If the only processes taking place is information being 
swapped by nearest neighbors, then as Stephen Wolfram 
proposes, there may be a universal rule set that governs nearly 
all the dynamics observed in the universe at all scales [4]. 

In this manuscript we will describe methods of exploiting 
the nearest neighbor physics of meso-scale phenomena for 
computation. As Yashihito [5], and many others have pointed 
out, as the device sizes shrink and the level of integration 
in microcircuits increases we will more closely approach the 
nanoscale. What was clearly articulated by Yashihito is that we 
should be able to use matter itself for our computations. We 
should be able to exploit the molecular dynamics and meso- 
scale physics for computations. Yashihito did not make explicit 
suggestions on how to undertake this task. Miller suggested a 
variety of physical systems that might be configured to carry 
out computation [6]. One of the suggestions was liquid crystal. 
This has recently been shown to be possible; Harding and 
Miller [7], [8] have demonstrated for the first time that liquid 
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Fig. 1 . Schematic of proposed computational system witli bulk matter 



crystal can be evolved to do analogue filtering. This will be 
discussed in section |lll] In the following sections we outline a 
number of approaches to exploit the physics of materials for 
computation. We describe numerical and/or empirical results 
for our suggestions. 

Figure [T] shows a schematic of the proposed technology. 
Basically, we will utilize a block of matter (solid, liquid, 
or gas) that allows us to change its properties/behavior by 
external forces. The external forces induce property/ changes, 
which we can consider to be a sort of "computer program.", 
so that there is a direct link between the external forces that 
we have control over and the induced changes in the block 
of matter. By measuring the behaviour of the altered block of 
matter, we can essentially submit input data to the sample and 
receive output data. In this way we have performed a type of 
computation. 

Of course we cannot directly program the molecular dy- 
namics and we do not have control of the molecules, at least 
not directly. The molecules will interact with their nearest 
neighbor and we can exploit these phenomena along with the 
state changes in regions of the block of matter, in order to 
perform computations. Figure [T| shows a nanoscale schematic. 

In the following sections we introduce the concept of 
programmable matter by first describing a well-known "pro- 
grammable surface" exploited by electrical engineers. Then we 
describe several examples of programmable matter. Our first 
example describes a technique that can be used to "program" 
liquid crystals to perform a computation in the form of signal 
processing. Following that we describe the use of a magnetic 
quantum dot array for associative memory. In the next example 
we introduce the idea of a programmable Fermi surface in cer- 
tain types of solid-state crystal lattices, in which we describe 
a prototype system and present some simulation results. In a 
final example we describe a technique to use acoustic pulses 
to modulate the density of bulk matter and thereby modulate 
the refractive index of the material. These refractive index 
modulations could be used for optical computing. 

We stress that some of the examples presented are in 
preliminary stages of investigation. Therefore, the paper is 
somewhat speculative but we believe the early results tend 
to support our speculations, and we believe the preliminary 
results will be of interest to a wider research community. 
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Fig. 2. Schematic of the architecture for an FPGA. Each cell can behave 
like any two-input, one-output Boolean logic function. 



II. Field Programmable Gate Arrays 



Electrical and computer engineers utilize a chip known as 
a Field Programmable Gate Array (FPGA). These chips are 
capable of being configured with software and can emulate 
many types of digital circuits. The architecture, schematically 
shown in Figure |2] is analogous to how we can exploit the 
physics of materials for computation, so a detailed discussion 
is important. 

The chip consists of an array of programmable logic cells 
that are connected to their nearest-neighbors. Each cell can 
exhibit any two-input one-output Boolean logic function. 

In addition the connections between the cells are pro- 
grammable so essentially any Boolean network can be con- 
figured within the space limitations of the chip. The chip 
essentially is a sea of gates - a programmable surface. 

Languages exist for programming these chips. In practice 
one develops a "circuit" which is then downloaded into static 
random access memory (SRAM). When a system utilizing this 
combination of SRAM and FPGA is powered up the SRAM 
sends configuration information to the FPGA and effectively 
"wires" up the logic cells. This download or wiring can take 
several milliseconds depending on the complexity of the circuit 
and the number of bits being downloaded. 

It is also possible to program the FPGAs as if they were 
black-boxes without any knowledge of the internal connec- 
tions. This knowledge-free heuristic programming approach is 
exactly the same technique we can use to exploit the physics of 
materials for computation. The technique is called evolutionary 
programming. 

Wolfram was one of the first to suggest using evolutionary 
algorithms for evolving real-world systems that may be too 
complex for humans to engineer [9]. 
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One of the first hardware implementations using evolu- 
tionary programming in an FPGA was described by [10] 
and [11]. They describe evolution of an analog filter in an 
FPGA. Significantly, they discovered that evolution designed 
circuits that were often irrational from an engineering design 
perspective. The circuits would exploit parasitic capacitances 
and inductances of the transistors. Thompson was the first 
to point out that evolution was exploiting the physics of the 
gate in order to achieve the desired fitness. Later work by 
him [12] was on evolving circuits (in simulation) with single 
electron transistors. These transistors are expected to operate at 
K, and so the simulations typically assume this temperature. 
Instead Thompson and Wasshuber assumed 340mK, which is 
too hot for the single electron transistor, and thermal noise 
becomes significant. What they found was that evolution was 
able to exploit the thermal noise to achieve the goal of building 
a NOR gate. Later simulation work by Thompson describes the 
evolution of built in self-test circuits for evaluating sequential 
and combinatorial circuits [13]. This was unique in that the 
self-test circuit actually utilized some of the same components 
used by the main circuit. Finally [14] describe evolution in 
FPGA circuits. They describe a tone discriminator very similar 
to [11]. Harding and Miller [7], [8] have recently shown 
that a similar technique can be used to exploit the molecular 
interactions in liquid crystal and they have evolved "circuits" 
that perform analog filtering(see section HIH i. 

In the following, we review some work of Huelsbergen 
et al concerning the utilisation of a genetic algorithm and 
an FPGA to evolve an oscillator circuit [15], [16]. In this 
work, they selected a region of the chip to be exploited, 
and then decided which pin would be the output where the 
oscillations could be observed. Recall that the cells and the 
connections can be configured. The genetic algorithm (GA) 
- a type of evolutionary algorithm - is easily described. In 
order to use GA to find an oscillator on this programmable 
surface (the sea of gates), the configuration information is 
coded as a long bit string of many thousands of bits. This 
string is called the chromosome or the genome. Different 
segments of the string (genes) code for specific configuration 
information in the FPGA. For example, one segment could 
configure the cells and another segment could configure the 
connections. In practice it is a little more complicated, in as 
much as different segments of the string would code for north 
connections, another segment of the string would code for the 
south connections, etc. The genome consisted of about 10,000 
bits. 

Since we know that the genome codes all the information 
to be sent to the FPGA and represents the entire configuration 
of the circuit, we then need some fitness metric to evaluate the 
said circuit. Huelsbergen et al. chose to monitor the pulses (or 
lack of pulses) on a selected output pin. But they needed to 
monitor the pin long enough to integrate over that time period, 
to determine if the circuit was oscillating, and determine the 
frequency of the oscillations. 

After selecting the genome - the bit string representing 




Fig. 3. FPGA - programmable surface - used by [15], [16] for their 
experiments in evolution in a physical realm. 



the circuit - and fitness - a metric of how well the circuit 
behaves - the next step in the GA consists of initializing 
a population of these strings and evaluating each of them. 
This consists of actually downloading each string in sequence 
and evaluating how well the circuit behaves with respect to 
the goal. The worst performing strings are discarded. The 
best strings are kept and segments between best performing 
strings are swapped with each other to create a new population 
of strings, which is then evaluated. The process repeats and 
eventually converges to strings/circuits that have either the 
desired behavior or is close to the desired behavior Figure 
[3] is a photo of the circuit they used for their experiments. 

There was some peculiarity in their results, which have 
impact on the use of evolution algorithms for configuring phys- 
ical media. Figure |4] is an example of the oscilloscope trace 
from one of the oscillators that was evolved. Firstly, notice 
that the waveform is a saw-tooth, not the expected square wave 
from a CMOS device. Secondly, notice the voltage levels. The 
chip can be configured at TTL or CMOS voltage levels. In this 
experiment it was configured at TTL and this is seen in the 
voltage levels for saw-tooth wave. The lowest voltage point 
is about 1.5 volts. This is the highest voltage that TTL logic 
will register logic low. The highest voltage, excluding spikes, 
is at 3.8 volts. This is the lowest that TTL logic will register 
high logic. This, along with the random voltage spikes above 
5.0 volts indicates the transistors are operating in the linear 
regime and relaxing by releasing random spikes. This bizarre 
behavior from CMOS logic is due to the fact that there were no 
constraints in the fitness function to force the logic to operate 
in the digital regime. Evolution does not know anything about 
the physical medium. It is blindly optimizing a fitness function 
given the physical medium. So the evolved circuit essentially 
exploits the physics of the transistors themselves without any 
specific instructions to do so. Furthermore, the evolved circuit 
is sensitive to the actual FPGA. If we replace the FPGA 
with another chip and download the evolved circuit it will 
behave slightly differently with respect to the frequency. This 
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Fig. 4. Oscilloscope trace from evolved oscillator circuit in CMOS hardware 
(Huelsbergen et al. 1998) [15]. 



candidate for evolution in materio would be liquid crystals 
[18]. Recently this suggestion has been supported by work of 
Harding and Miller [8], who showed that it is relatively easy to 
configure (using computer controlled evolution) liquid crystals 
to perform various forms of computation. 

For this work we have chosen to use a genetic algorithm, 
as there is an already well established field of evolvable 
hardware from which we can draw experience. Other search 
methods have been suggested for material systems, and will 
probably work to the same level. In previous work, Toffoli 
argues that simulated annealing may provide a more suit- 
able programming technique for programmable materials [19]. 
This technique shares many similarities with evolutionary 
algorithms, however, simulating annealing can perform less 
efficiently as it does not use a population based approach and 
therefore can more easily become trapped in local attractors. 



is further evidence that evolution is exploiting the physics of 
the transistors. When we replace the transistors with a different 
physical set of transistors (i.e. replace the silicon chip) we get 
different behavior. 

In summary, we have shown that we can use evolution to 
design computational circuits with a "sea of gates." Evolution 
will exploit the substrate at the lowest level possible to 
accompUsh the desired task. Evolution will actually exploit the 
physics of the gates and the variation in the silicon crystals to 
achieve its goals. There are several implications for exploiting 
the physics of matter Firstly, we do not need to pre-compute 
anything. Secondly, by definition, since we are manipulating 
physical matter with the evolutionary algorithm if a solution 
can be found then we will find it and it will be physically 
realizable. Thirdly, evolution will cobble together any thing 
to get the problem solved. It will not necessarily "build" the 
"computational system" we expect, because it will exploit the 
physics of the medium in ways we may not even imagine. 



III. Programming a Liquid Crystal Substrate 
A. Intmduction 

As shown above in Figure |4] and discussed by Thompson et 
al. [17] and Huelsbergen et al [15], [16] it is now obvious that 
evolutionary algorithms exploited subtle physical properties of 
the FPGA and associated circuits in order to solve a particular 
problem. It is not fully understood what properties of the 
FPGA substrate were being exploited. This lack of knowledge 
of how the system works prevents engineering the design of 
systems to exploit the complex physical characteristics. We 
argue that the lesson that should be drawn from the work 
of Thompson is that evolution may be used to exploit the 
properties of a wider range of materials than silicon. This 
further suggests an exploration through artificial evolution. We 
refer to this as "evolution in materio." Miller suggested a good 



B. Liquid Crystal 

Liquid crystal (LC) is commonly defined as a substance 
that can exist in a mesomorphic state [20] [21]. Mesomorphic 
states have a degree of molecular order that lies between that 
of a solid crystal (long-range positional and orientational) and 
a liquid, glass or amorphous solid (no long-range order). In 
LC there is long-range orientational order but no long-range 
positional order 

Aromatic LC is often called a benzene derivative. There 
is also heterocyclic LC where one or more of the benzene 
rings are replaced with pyridine, pyrimidine or other similar 
group. LC can also have metallic atoms (as a terminal group) 
in which case they are called organometallic compounds. 
Chemical stability is strongly influenced by the linkage group. 
Compounds where the aromatic rings are directly linked are 
extremely stable. LC tends to be transparent in the visible and 
near infrared and quite absorptive in UV. 
There are three distinct types of LC: lyotropic, polymeric and 
thermotropic. Thermotropic LC (TLC) is the most common 
form and is widely used. TLC exhibit various liquid crystalline 
phases as a function of temperature. They can be depicted as 
rod-like molecules and interact with each other in distinctive 
ordered structures. TLC exists in three main forms: nematic, 
cholesteric and smectic. In nematic LC, the molecules are 
positionally arranged at random, but all share a common 
alignment axis. Cholesteric LC (or chiral nematic) is like 
nematic however they have a chiral orientation. In smectic LC, 
there is typically a layered, positionally disordered structure. In 
type A the molecules are oriented in alignment with the natural 
physical axes (i.e normal to the glass container, depicted by 
the arrow), however in type C, the common molecular axes of 
orientation is at an angle to the container. 

There is a vast range of different types of liquid crystal. LC 
of different types can be mixed. LC can be doped (as in Dye- 
Doped LC) to alter their light absorption characteristics. Dye- 
Doped LC film has been made that is optically addressable and 
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can undergo very large changes in refractive index [22]. There 
are Polymer-Dispersed Liquid Crystals these can have tailored 
electrically controlled light refractive properties. Another in- 
teresting form of LC being actively investigated is Discotic 
LC. These have the form of disordered stacks ( 1 -dimensional 
fluids) of disc-shaped molecules on a two dimensional lattice. 
Although discotic LC is an electrical insulator, it can be 
made to conduct by doping with oxidants [23]. LC is widely 
known as useful in electronic displays, however, there are 
in fact, many non-display applications too. There are many 
applications of LC to electrically controlled light modulation: 
phase modulation, optical correlation, optical interconnects 
and switches, wavelength filters and optical neural networks. 
In the latter case, LC is used to encode the weights in a neural 
network [24]. 




Fig. 5. Equivalent circuit for LC 

Figure |5] shows the equivalent electrical circuit for liquid 
crystal between two electrodes when an AC voltage is applied. 
The distributed resistors, R, are produced by the electrodes. 
The capacitance, C, and the conductance, G, are produced by 
the liquid crystal layer [25]. 




Fig. 7. Pliotograpli of the LCEM prototype circuit. 



addressing pixels so that images can be displayed. The driver 
circuit has a large number of outputs that connect to the wires 
on the matrix display. When displaying an image, appropriate 
connections are held high, at a fixed voltage - the outputs are 
typically either fully on or fully off. 

Such a driver circuit was unsuitable for the task of intrinsic 
evolution. There is a need to be able to apply both control 
signals and incident signals to the display, and also record 
the response from a particular connector. Evolution should 
be allowed to determine the correct voltages to apply, and 
may choose to apply several different values. The evolutionary 
algorithm should also be able to select suitable positions 
to apply and record values. A standard driver circuit would 
be unable to do this satisfactorily. Hence a variation of the 
evolvable motherboard was developed in order to meet these 
requirements. 



C. An Evolvable Motherboard with a FPMA 

We have been experimenting with a hardware system that 
enables programmability and reconfiguration through a semi- 
conductor cross point array switch. We refer to this system 
as a evolvable motherboard (EM) [26]. When the EM is 
connected to a material substrate (e.g. liquid crystal display) 
or discrete electronic components [26], [27] we refer to the 
system as a field programmable matter array (FPMA). The EM 
is connected to a PC that is used to control the evolutionary 
processes. The EM also has digital and analog I/O that can be 
accessed for test and recording of the response of the material 
under evolution. 

In the experiments presented here, a standard liquid crystal 
display with twisted nematic liquid crystals was used as the 
medium for evolution. It is assumed that the electrodes are 
indium tin oxide. Typically, such a display would be connected 
to a driver circuit. The driver circuit has a configuration bus on 
which commands can be given for writing text or individually 



The Liquid Crystal Evolvable Motherboard (LCEM) is a cir- 
cuit that uses four cross-switch matrix devices to dynamically 
configure circuits connecting to the Uquid crystal. The switches 
are used to wire the 64 connections on the LCD to one of 
8 external connections. The external connections are: input 
voltages, grounding, signals and connections to measurement 
devices. Each of the external connectors can be wired to any 
of the connections to the LCD. 

The external connections of the LCEM are connected to 
the computers analogue inputs and outputs. One connection 
was assigned for the incident signal, one for measurement and 
the other for fixed voltages. The value of the fixed voltages 
is determined by the evolutionary algorithm, but is constant 
throughout each evaluation. Each of the external connectors 
can be wired to any of the connections in the LCD (see figures 
|6]and[8]l. 

In these experiments the liquid crystal glass sandwich was 
removed from the display controller it was originally mounted 
on, and placed on the LCEM. The display has a large number 
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characteristics of the LCD are known. This raises the possibil- 
ity that a configuration may be appHed that would damage the 
device. The wires inside the LCD are made of an extremely 
thin material that could easily be burnt out if too much current 
flows through them. To guard against this, each connection 
to the LCD is made through a 4.7Kohm resistor in order 
to provide protection against short circuits and to help limit 
the current in the LCD. The current supplied to the LCD is 
Umited to 100mA. The software controlling the evolution is 
also responsible for avoiding configurations that may endanger 
the device (such as short circuits). 

It is important to note that other than the control circuitry 
for the switch arrays there are no other active components on 
the motherboard - only analog switches, smoothing capacitors, 
resistors and the LCD are present. 



D. Genetic Representation 



The genetic representation for each individual is made of 
Fig 8 Scliematic of LCEM two parts. The first part specifies the connectivity; the second 

part determines the configuration voltages applied to the the 
LCD. 



of connections (in excess of 200), however because of PCB 
manufacturing constraints we are limited in the size of con- 
nection we can make, and hence the number of connections. 
The LCD is therefore roughly positioned over the pads on the 
PCB, with many of the PCB pads touching more than 1 of 
the connectors on the LCD. This means that we are applying 
configuration voltages to several areas of LC at the same time. 

Unfortunately neither the internal structure nor the electrical 



Each connector on the LCD can be connected to one of 
the eight external connectors or left to float, figure [8] Each of 
the connectors is represented by a number from to 7 and 
no connection is represented by 8. Hence the genotype for 
connectivity is a string of 64 integers in the range to 8. 

The remainder of the genotype specifies the voltages applied 
to the pins on the external connector that are not used for 
signal injection / monitoring. On the LCEM there five such 
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configurable connectors - the other three are connected to 
ground, the incident signal and to the data recorder. The 
voltage is represented as a 16-bit integer, the 65536 possible 
values map to the voltage levels output from -lOV to +10V. 
The second section of the genotype is therefore represented as 
a string of five 16bit integers. 

A mutation is defined as randomly taking an element 
in one part of the genotype and setting it to a randomly 
selected new value. Constraints are enforced to prevent illegal 
configurations. 

We chose not to use genetic recombination as the constraints 
imposed on this representation would make it difficult to 
implement and would require many arbitrary decisions to be 
made on suitable repair techniques. For example, it is unclear 
what strategy should be used to fix a genotype where there are 
multiple outputs and only one is allowed. For this reason, the 
evolutionary algorithm used here has no crossover operator. 

In all the following experiments, a population of 40 individ- 
uals was used. The mutation rate was set to 5 mutations per 
individual. Elitism was used, with 5 individuals selected from 
the population going through to the next generation. Selection 
was performed using tournament selection based on a sample 
of 5 individuals. 

Evolutionary runs were limited to 100 generations. With 
each generation taking approximately 60s to evaluate. 

E. Evolution Of A Tone Discriminator 

A tone discriminator is a device, which when presented 
with one of two input signals returns a different response 
for the each signal. In [17], on which this experiment is 
loosely based, the FPGA under investigation was asked to 
differentiate a IkHz square wave from a lOkHz square wave, 
giving a low output for one and a high output for the other. 
In this experiment we have arbitrarily chosen two frequencies 
of lOOHz and 5kHz. Each signal is a square wave, oscillating 
between OV and 5V, with equal timing given to the low and 
high states. The tones were presented in 250ms bursts with no 
gap between the tones. The goal was to evolve a device that 
would output a low value ((O.IV) at low frequencies, and high 
()0.1V) at the higher frequency. The fitness was calculated as 
the percentage of samples made where the output was in the 
correct state for a given input frequency. 

Let S be the vector containing the input sample. Let L be 
the length of S". O is a vector containing the output frequency 
at a given time. The output frequency can be either HIGH or 
LOW. The jth element of the set is S [j] . t is the threshold 
for a low response, i.e. (O.IV. 

( 1 if S[i\<t and O [i] = HIGH 
x{i) = < 1 if S [{\>t and O [i] = LOW 
I otherwise 
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Fig. 9. Tone discriminator response. Dai'k areas indicate 5kHz input, liglit 
lOOHz 

fitness = 2^0^ 

It is important to note that samples taken do not correspond 
to time as samples are taken on an interrupt, and frequency of 
sampling may be affected by other processes running on the 
computer. 

Not all attempts at evolving the discriminator under these 
conditions were successful, however we did manage to evolve 
a discriminator with the response shown in Figure |9] Although 
the output was not stable, there is a clear difference between 
the behaviour at low and high frequencies. At high frequencies, 
a high output was obtained for the majority over time and for 
low frequencies a low output was obtained. 

We assume the behaviour stems from capacitative effects 
originating inside the LCD, and that the system is acting 
as a form of R-C network. The crosspoint switches are 
unlikely to be involved as they are designed for high frequency 
audio/video signals. The feed-through capacitance at IMhz is 
0.2pF and the switch I/O capacitance is 20pf. This would seem 
too small to have any filtering effect on these relatively low 
frequencies. 

An interesting observation is that if a configuration is 
reloaded into the LCEM it fails to work, however if the 
population containing that solution is allowed to evolve (for 
another 2 to 3 generations) the behaviour returns. The cause of 
this appears to be a lack of stability in the system, however it 
is unclear if this is caused by the liquid crystal itself or some 
other component in the system. 



F. Summary 

The work described is the first example of the use of com- 
puter controlled evolution to build computational processors in 
liquid crystal. Although liquid crystal appears to be suitable 
for use as an evolutionary medium, there are many unanswered 
questions. More work is required to prove that the LC is 
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responsible for the observed results and to attempt to discover 
how the LC is being exploited. It may be that the LC layer 
is being used as some form of configurable continuous RC 
network - similar to that shown in figure |5] but with much 
richer properties. 

IV. Magnetic Quantum Dots 

There are significant parallels between neural networks, 
chaotic computing, optical computing, molecular computing, 
and of course, quantum computing [28]-[33]. All of these 
have in common the fact that they are "unconventional models 
of computation." Two of the most interesting computational 
paradigms that join several of the other unconventional ap- 
proaches are - quantum dot cellular automata and quantum 
spin computing. 

Benioff first suggested quantum mechanical computation as 
Turing machine automata on Hamiltonian lattices [34], [35]. 
Albert suggested a similar concept of quantum mechanical 
automata [36]. These really significant works no doubt in- 
spired the suggestion of cellular automata for nanometer-scale 
computing [37], and inspired Porod et al. for their work on 
quantum-dot cellular automata [38]. (Cellular automata will 
be introduced shortly.) 

Other quantum computing work has drawn parallels with 
brains. The most notable in this effort is Penrose [39], [40] 
and the more recent is Satinover [41]. Though these works 
are largely in the area of philosophical research; it has had 
some impact on more serious physics research such as Kak 
[42] who wrote about state generators and complex neural 
memories in neural networks. This early Kak paper resulted in 
his later work on comparisons with brain dynamics and neural 
networks [43], which is typical in the field. Finally, Ventura 
implemented competitive learning in a quantum computing 
neural network [44] ; Behrman described a spatial quantum 
neural computer [45]; and Zak described quantum analog 
computing [46]. Though most of this work is very theoretical 
it an important contribution. Of a more experimental nature, 
others have described methods for building quantum networks 
[47]. Much of that approach suggests quantum spin systems 
using quantum dots. 

One of the earliest suggestions of quantum dot spin systems 
was described by Smith [48]. Their suggestion was to build 
arrays of GaAs-based quantum dots. A review of related 
experimental work is given by Chakraborty [49]. Spin tran- 
sistors and spin electronics are finally attracting a great deal 
of attention ( [50]-[52]). 

We propose an experimental prototype and describe simu- 
lation results, which are a follow up from some interesting 
observations by Kirczenow et al [53] on the dynamics of 
quantum hall phenomena in arrays of magnetic quantum dots, 
and some suggestions of [54] on quantum computing using 
dipole-dipole block systems. The work by Kirczenow et al. 




Fig. 10. Ising spin array for quantum computing 

describes two-dimensional arrays of magnetic dots with larger 
dots on the edges that are used to couple into the array. They 
show that if the dot is larger in diameter than in height it 
will take on only two spin states. In this respect these results 
are similar to [55] and [56]. In the following we describe 
some simulations we have done to show that it is possible to 
use a genetic algorithm to program a magnetic quantum dot 
computer for use in pattern recognition. 



A. Ising Spins 

In this section, the physics of spin glasses or Ising models, 
follows the presentation by [57]. [58] is a good introduction 
to the statistical physics of spin glasses. The subject of Ising 
spins has recently been re-introduced by [59]. 

Consider a square grid of lattice points. At each point we 
will place a tiny magnet (a spin) and we will represent this 
by spin up or spin down. These spins are represented by two 
colors, black and white, in Figure[TO] Now let the spins flip so 
the entire system can relax to the lowest energy configuration. 
This can be done at each temperature and it turns out there is a 
minimum energy for a given grid of spins at any temperature. 
At high temperatures the spin arrangement can be highly 
random with little or no correlations. At lower temperature, we 
will find that the spins are more spatially correlated. That, in 
summary, is the entire idea behind Ising models for magnetic 
spins (a.k.a. spin glasses). 

In the above, the energy is a function of the spin, E [{si}], 
for each spin in the system. So the total energy is given by 

where ei{si) is the energy of an individual spin that is not 
influenced by another neighboring spin. Since our system is 
binary we can write this equation as 



E[{s^}] = iE.h(l)-e,(-l)]s,+ 
K(l) + e,(-l)] 

^Eo-E^h^S^ (1) 
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where some obvious terms have been collected together as 
the "ground-state" energy, Eq . The quantity hi is the energy 
due to orientation of the spins, and in a magnetic system they 
represent the external magnetic field. This external magnetic 
field will prove to be very useful in programming an Ising 
spin computer 

In this simple system we can write the probability of a 
particular configuration of spins by the following equation; 



J~J . exp{hiSi) 



We let a "partition function" be given by: 

^ = E.ea;p(-i?N)/T 

then we can find the minimum energy by computing the 
partial derivative of the natural logarithm of Z with respect 
to the inverse of temperature (3. 



HZ) = \n{J{exp{|3h^) + exp{-[3h^)^ 

= ^ ln[exp{l3hi) + exp{—l3hi)] (2) 




Fig. 11. 



eisrgy^ 

Magnetic domain "transfer function" 



We now expand the model to allow the spins to interact with 
their neighbors. The interaction energy, J, between spins may 
come about through several mechanisms, but the end result 
will be that the spins will take on either a direction J > or 
J < 0. In a similar way to the derivation of the energy relation 
of Eq [T] we can write an energy relation including the spins 
and interaction energy, J and get the following equation: 



From this we can write the minimum energy as 



Sln(Z) 

_ hi{exp{phi) - exp{-(3hi)) 
~ ^ [expifihi) + exp{-ph{)) 

= ^h, tanh(/3/i,j) (3) 

■i 

Figure [TT] is a plot of this function. 

This is a very interesting representation for the magnetic 
field, because it parallels computational systems known as 
neural networks (cf. [60] and [61]). We can write the mag- 
netization TO as 



m = (sj) = tanh(/3ft,i) (4) 

This equation says the hyperbolic tangent of spin energy 
divided by the temperature gives the magnetization. For neural 
networks the output of a neuron is given by the hyperbolic 
tangent of the sum of the inputs to that neuron. As is 
well known the hyperbolic tangent function is graphically 
represented by a function exactly like that shown in Figure 

m 



J T,ij Si+lj + Sij S^j + l ) (5) 



This equation includes more dimensions because now we 
are including the neighbor interactions. In the previous equa- 
tion we did not include neighbors so the current equation 
appears to be more complex. 

Briefly, here is how we can use Ising spins for computation. 
We wiU view a magnetic quantum dot, Ising array as a dynam- 
ical system that is carrying out some specific computational 
task. For example, by associating some input state with a 
reproducible output state we can build an associative memory. 

We propose using an array of magnetic quantum dots that 
have locally programmable regions of the dot array. Various 
regions can be programmed with magnetic probes to induce 
specific magnetic spins in the dots. These dots can then interact 
with their nearest neighbor as a complex dynamical system. 
By using an array of magnetic probes on the edges as input 
and output channels, we can exploit the entire landscape of 
magnetic quantum dots for analog computation. Furthermore, 
the system is fault tolerant. If segments of the array are 
damaged, or nonfunctional, the array is dynamic and can adapt 
to these defects. 
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t=n+1 



the computation will be adjusted by controlled perturbations 
of the quantum dots. These perturbations are equivalent to 
programming the system. And the final spatial configuration 
is the computed answer to some mathematical mapping. 

So summarizing the dynamics of the Ising computer we can 
rewrite the relation for energy of the magnetic interaction as 



Fig. 12. Spin interaction map 



-{Br 



JijSiSj) 



(6) 



B. Simulation of Magnetic Dot Computing 

A cellular automaton is an array of small computational 
elements called finite state machines. Each cell in the array 
is capable of receiving input from its nearest neighbors only. 
At each time increment (time is discrete), the input from the 
neighbors determines the new state for the cell. For example, 
consider the small array in Figure [12] and the cell marked C 
in that array. At time t ~ n we see that cell C is white and 
is surrounded by two black cells. At time t = n + 1 we see 
the center cell is now changed to a black state. This change 
is due to dynamic interactions between the four neighbors 
and the center cell. Every cell in the entire array is updated 
and the updating can be synchronous or asynchronous. If 
synchronous, than the next state for each cell is computed 
before any changes take place in the array. At that time all 
the cells undergo the computed state change. If the dynamics 
are asynchronous then the computations are done for a cell 
selected at random and the cell undergoes the required change. 
This in essence, is a cellular automaton. 

Both electronic quantum dot arrays ( [38]) and a magnetic 
quantum dot arrays are analogues to cellular automata. In 
each case, the computations carried out at each "dot" is done 
by nearest neighbor interactions. The electronic quantum dot 
array can be used for computational purposes by electronic 
interaction between nearest neighbor quantum wells. The 
magnetic quantum dot array can be used for computational 
purposes by nearest neighbor interaction of the magnetic 
domains. All of these systems have been modeled using 
cellular automata. 

As pointed out in [62] and [63], the primary problem associ- 
ated with using cellular automata or dot arrays for computation 
is the fact that the flow of information is bidirectional. There 
is no isolation between the inputs and the outputs. Here we 
present an approach suggested in [63], but we include enough 
detail to make the suggested methods practical, and we present 
simulations to demonstrate feasibility. Though we expect the 
method will also work with regular electronic quantum dots, 
cellular automata, and molecular oscillator-based computing. 
The key to utilizing these systems for computational purposes, 
is to recognize that the dynamics of the system are sensitive 
to initial conditions. This is equivalent to setting up the input 
data for a computation. The evolution of the dynamics of 



That is, the energy at each site is given by the product of the 
neighboring spins and the interaction energy J. Here we have 
also included an external magnetic field B. The probability of 
a given spin state at some site is given by 



exp 



- (E, J^jSiSj/T 



(7) 



This has the same functional form as the plot shown in 
Figure [TT| The figure is a plot of the probability at several 
temperatures. At very low temperatures, the probability will 
tend to take on either zero or one. At warmer temperatures 
the probability will tend more-and-more toward P = 0.5. 
These two equations (Eq II6I7II ) are the only ones needed for 
a simulation of the Ising Computer. 

Typically in running simulations of Ising systems one would 
use periodic boundary conditions. Wrapping the edges will 
accomplish this, so that the top and bottom are connected, 
and the left and right edges are connected. Since we want 
to use the system for computation we will use the edges for 
input and output, so we will not use the periodic boundary 
condition. Figure [13] shows the spins for a 20X20 system at 
various snap-shots in time. We see from Figure [13] that as 
the dynamics progresses in time the spins become more-and- 
more locally correlated. The effects of temperature are shown 
in Figure [14] As indicated in Equation [7] and Figure [14] the 
temperature affects the dynamics. 

To demonstrate the programmability of the system we used 
three edges as the input and one edge as the output. We 
used a genetic algorithm (cf. [64]) with three, 2-dimensional 
chromosomes per member in a population of Ising systems. 
The three chromosomes represent the initial configuration, 
/[x][y], the coupling magnetization between sites, J[x][y], 
and an added magnetic field, i3[x][y], to increase the dynamic 
programming ability of the system. The relevant dynamics is 
given by the probability 



P. = 



1 + exp 



(8) 



This form of the equation hardly differs much from Equa- 
tion [7] The only difference is the inclusion of the external 
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Fig. 13. Time dynamics of tlie magnetic quantum spins. 





Fig. 14. Temperature study of magnetic spin quantum dot computer. 
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Fig. 15. Learning curve for pattern matching problem using magnetic 
quantum dot computer 



magnetic field B. 

Figure [15] shows the results from a simulation. The system 
works as follows. The genetic algorithm manipulates the 
elements of the three matrices, /, ,/, B such that a particular 
input pattern can be mapped to a particular output pattern 
with a low Hamming distance between the actual output and 



the target-output pattern. In a sense, the system acts like a 
feedforward neural network. The information on the input 
edges is dynamically interacting with their nearest neighbors. 
These neighbors in turn are interacting deeper into the array. 
So the information from the input edges is essentially fed into 
the array. The array relaxes to some configuration with the 
output edge giving the final computation. 

In the initial experiments all three matrices (/, J, B) were 
initialized by the genetic algorithm and the system was allowed 
to settle to the answer, some pre-selected pattern on the output. 
It should be possible to also use more complex arrays of 
J{t) or B{t) and to let these represent a dynamic program. 
For example, with an array size of 20X20 one would have 
an active dynamic array of 19X19, assuming all edges are 
used for input/output, and then the arrays J and B would be 
configured as follows. J[a;][?/][t] and i?[a;][y][t] where each 
x,y array also has a time index from t = \ lo t — 19. This is 
easily computed by the fact that the signal can only travel one 
spatial site per time step (called the speed of light). Of course 
it should be possible to use any t index depending on the 
problem. For example, some problems may require only t — n 
updates of the B and J matrices, while others may require 
more. It could also be possible to run the dynamical system 
until n time units after the B, J updates to the program. In 
this case the free running time after set-up, would be part of 
the programming of the system. 
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C. Prototype System 



Programming array 



In the above simulation we made several assumptions, 
primarily concerning size of the magnetic dots and how the 
dots would interact with each other First, we assumed that the 
dots would be larger in diameter than in height. This important 
assumption will result in the dots taking only one of two 
states (spin up and spin down). Shinjo et al. first observed 
this experimentally (2000) [65]. Their observations indicate 
that dots on the order of microns in size (1-2 microns) will 
still exhibit only two states. Dots this large suggest that we 
can make a prototype easily. 




inpu 



oi.tpurs 



Coils around 
inputs pin's 

Fig. 16. Scliematic of a prototype for magnetic quantum dot computing. 



We know that the magnetic quantum dots would not really 
be able to implement quantum computation. The magnetic dots 
are too large to be called quantum dots. Classical electrostatics 
will dictate the overall dynamics. Cowburn et al. describe room 
temperature quantum cellular automata made with magnetic 
quantum dots [66]. When a cluster of dots takes on the same i 
spin state it acts essentially as a single giant classical spin 
system. Their quantum dots were about 110 nm in diameter 
and 10 nm thick. I ^.^ ' 

As demonstrated in [65] and [66], small magnetic fields 
are needed (e.g. 300 Oe) for programming these dots. Both 
sets of investigators use small magnetic "pins" to inject fields 
into specific locations. Shinjo et al. uses a magnetic probe to 
monitor the spin state and Cowburn and Welland use reflected 
light to access the longitudinal Kerr effect to observe the spin 
state. 

Figure [16] shows a schematic of a prototype system that 
could easily be constructed. Some of the edges would be inputs 
and another edge would be the outputs. The edge probes are 
held fixed in position and the inputs are held at a positive or 
negative voltage so as to maintain the input state to the array. 
The output probes continually monitor the output dots. In the 
center of the array is a small cluster of magnetic probes that pjg 17 
raster back-and-forth to program the initial state. 




r!i.h 



Schematic of a programmable resistor aiTay. 



Using this technique to raster a cluster of probes, in which 
the probe size is potentially larger than the magnetic quantum 
dots, we can still perform useful computations. In our simula- 
tion we assumed that we could not necessarily access each dot, 
but rather a small cluster of dots with each magnetic probe. 
In other words, the probe size was larger than the dot size. 

In [66] and [67] a technique is described that uses RAM 
cells to program the individual dots. They calculate that a 
microprocessor based on magnetic quantum dots would draw 
only one Watt of power Of course, electronic control of the 
dots is also possible, for example spin transistors [67] and 
[68]. 



V. Programmable Fermi Surfaces 



Consider a programmable resistor array. It can be used to 
perform the mathematical operation of vector-matrix multi- 
plication. This operation is at the heart of target recogni- 
tion, pattern recognition, data compression, process control, 
associative memory, content addressable memory, and many 
other computations. In fact, vector-matrix multiplication is the 
primary mathematical operation that digital signal processing 
(DSP) chips undertake. 

To exploit a programmable resistor array for associative 
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Fig. 18. Crystal structure of AfoOg (side view, top view) 



memory (a type of neural network), for example, we recall 
that an associative memory is a mapping between two different 
vectors. When input nodes in the circuit are activated the cir- 
cuit will automatically perform a vector-matrix multiplication 
and the output nodes will present the nearest mapping relations 
to the input vector The matrix of connections can be pre- 
computed so that the system will respond to several (or many, 
depending on the system size) different input stimulates. 

The same system can be used for pattern recognition. The 
target pattern will be coded as a vector Then compute the cross 
product of that vector and its transpose. This will generate 
a matrix. This matrix is the resistance values for the circuit 
shown in Figure [17] - the programmable resistor array. When 
vectors similar to the expected target are presented to the 
inputs the circuit will compute a vector matrix multiplication 
and complete the recognition of the target vector (or image). 
All these computations are well known and fall into the 
category of neural networks ( [69]; [70]). These types of 
computations are also exactly the type that can be done with 
a programmable crystal lattice. 

Certain materials crystallize into stacks of 2-dimensional 
layers (Figure [TSl l. One example of this class of materials is 
molybdenum trioxide, an insulating material. It is possible to 
electrochemically reduce the molybdenum while it is in the 
crystal lattice and intercalate ionic species, such as A^a+, 
or Cs+ into the lattice (Figure [T9]l. The material will then 
become a semiconductor and fast-ion conductor. 

Kq ^MoOs is an excellent example of a fast-ion conductor 
Pure M0O3 is white to yellow in color. The potassium- 
doped crystal is dark red to purple in color; this darker color 
indicates that the electronic conductivity has increased. The 
doped crystal is a mixed ionic and electronic conductor. It is 
electronically a semiconductor and a potassium ion conductor 
The potassium ions disrupt the electron density (the Fermi 
surface) of the molybdenum atoms. Of course the electron 
density is an image of the charge density. 

In an ac field the ions will drift back-and-forth at the 
frequency of the applied field (Figure \19[ . So in a frequency 
and amplitude modulated ac field, the ions can be moved 
to within a few nanometers of each other to form clusters 
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Fig. 19. Effects of no field, dc field and modulated ac field on the ions in 
the lattice. 



and then "pinned" into place with a higher frequency field. 
Naturally these clusters of pinned ions would prefer to move 
away from each other and settle into some equilibrium state 
commensurate with the underlaying molybdenum oxide lattice. 
These metastable clusters of ions disturb the underlaying 
Fermi surface of the molybdenum ions, and generate waves of 
disturbance with a correlation length on the order of 100s of 
microns. These are known as charge density waves and they 
can be measure by their effects on ac and dc fields. Basically, 
by moving these ions in an ac field and pinning them we are 
in effect programming Fermi surfaces. 

A. Simulations 

In the absence of an electric field the ions will repel each 
other to a maximum distance and form a superlattice that 
may be commensurate with the under-laying primary lattice. 
In the presence of electric fields the ions will move and the 
superlattice will now become incommensurate. The lattice 
vibrations in this incommensurate region of the crystal will 
send out quasi-periodic or even chaotic waves ( [71]; [72]). 
It is believed that each of these ionic centers forms a phase 
vortex ring ( [73]). These vortex rings will expand until either 
a defect or another vortex ring is encountered. They can have a 
diameter of up to 100s of microns depending on the number- 
density of vortex rings ( [74]; [72]). The rings behave like 
quasi-periodic oscillators, known as sine-circle maps ( [75]; 
[76]; [77]; [78]; [72]). Obviously, a molecular scale array 
of them is a type of cellular automata known as a coupled 
map lattice. The number-density, which can be controlled by 
the stoichiometry or by the electric field modulation, will 
essentially dictate the size of the rings. We can exploit these 
phenomena for computation. 

We can "program" the vortex ring array by moving the 
ions (the vortex ring generating entities) with modulated ac 
fields. The ions can then be "pinned" into place with high 
frequency fields - about 10 MHz ( [79], [80]. The rings can 
be "connected" with external fields ( [81]; [72]; [82]; [73]; 
[83], [84]). At a specific threshold the rings will produce a 
sliding charge density wave. The actual "programming" can 
be done with a genetic algorithm searching for a "circuit" that 
computes a predefined mapping relation (exactly as was done 
with a FPGA search for oscillators). The best analogy is to 
think of this as nanoscale programmable resistor array - or 
programmable crystal lattices. 

Several interesting models of this phenomenon have been 
discussed in the literature. The simplest model is based on 
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Fig. 20. Myers and Sethna model 



Hamiltonian dynamics of particles trapped in wells ( [85]). The 
model starts by assuming a Hamiltonian lattice to describe the 
effects of an applied electric field. The charge density wave 
elasticity and the pinning disorder are given by 



^E^cos[27^(a;,-A)] (9) 

where Pi is a random number between and 1 selected for 
each site, i, but not at each time iteration. F is the applied 
field; V is the pinning potential, Jis the coupling strength (like 
a spring constant); and Xi is the phase at site i. The dynamics 
evolves according to 

The equation has solutions at zero when F < V . The 
solutions are shown graphically in Figure |20l The figure shows 
the phase and the phase velocity. The dynamics suggests a 
stick-slip behavior. The charge density wave gets stuck at some 
pinning center, and when the field increases, it will rapidly 
unstuck or slip and then get stuck at another point. 

The Myers and Sethna model is very interesting and has 
some similarity to the models of [75]; [76]; [77]; [78]; [72] 
who discuss sine-circle maps. Specifically Azbel and Bak 
suggest a model of the following type. 

Xn+i = Xn + ^ + {h/2t:) sin(27rx„) 

Azbel and Bak did not assemble array models of this 
equation where each cell in the array was a separate vortex 
interacting with their neighbors. In order to do that we will 
draw on the work of Kanek, who describes a simplified version 
of this equation in a cellular array known as a coupled map 
lattice (like a real-number cellular automata) [86]. We will 
make further assumptions based on the work of [87] who 
worked with coupled map lattices, and had the objective of 
demonstrating computation. In brief, we will combine the 




K 

Fig. 21. Bifurcation diagram for circle map. 



equation of Kaneko with the algorithmic modifications (to 
be described subsequently) of Sinha and Ditto and model 
computation with these sine-circle maps. 

The Kaneko modification of the sine-circle map (hereinafter 
called circle map) is given by 

Xn+i ~ x„ + (K/27r) sin(27rx„) 

We further simplified this equation as follows: 

Xn+l = Xn + K sill{NXn) 

The "bifurcation diagram" shown in Figure |2T| outlines the 
dynamics of this equation. Here we let TV = 1. From the 
bifurcation diagram it is clear a whole range of dynamics is 
possible. Zero is the attracting point from K = —2 to K = 0. 
A phase shift of occurs at K = and remain at that level 
untill K ^2. 

If each vortex site in the crystal lattice behaves like this 
model we can simulate an array of them and couple them 
together. This is called a coupled map lattice. To couple them 
together we will use an external threshold, T„ at each site. 
Connecting a linear array of these to a DC threshold we obtain 
the results shown in Figure l22l (see [87] for algorithm details). 

In this simulation the cells were operating with identical 
conditions and were in the chaotic regime. Figure l22h shows 
the chaotic signal and a threshold line. Whatever signal is 
above threshold will carry over into the next cell in the linear 
array (in a feed forward manor) and will add to the equation. 
Doing this, and monitoring the pulses at the last node in a 
chain of 100 cells, we find that periodic pulses develop after a 
short time interval. The figure on the right (Figure l22b) shows 
the periodicity. This is exactly one of the defining behaviors for 
charge density wave materials - periodic pulses are observed 
in a DC field ( [80]; [79]). 
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Fig. 23. Two-dimensional lattice simulation of charge density waves using 
circle-map equation and threshold dynamics of Sinha and Ditto (1998) [87] 
to represent an applied dc field 



With a two-dimensional array and two equal but perpendic- 
ular dc fields we expect to see traveling pulses above threshold. 
Figure |23] shows this phenomenon. The figure shows a two- 
dimensional plot of the above threshold values coming out of 
the oscillators. Periodicity in two-dimensions is clearly seen. 
It represents the "stick-slip" phenomena and gives rise to the 
sliding charge density wave. 

Sinha [87]-[90] demonstrates, in simulation, that chains of 
logistic functions can be used for Boolean logic and numerical 
calculations. For example, it is possible to build an adding 
machine by setting individual thresholds at each machine, and 
that the emitted signal at the end of the chain is equal to the 
addition of the thresholds. The mapping relation encoded the 
numbers and threshold is given in Figure |24] for both chains 
of logistic and circle map functions. Further, as to be expected 
arrays of these oscillators can be configured for constructive 
and destructive interference of pulses to produce Boolean logic 



functions. 



B. Prototype System 

Single crystals and epitaxy films will exhibit the phe- 
nomenon of charge density waves that can be exploited for 
analog computation by pinning the ions in specific locations. 
These materials were first investigated for solid-state batteries 
as an electrolyte, because they are ionic conductors. Obviously 
prototype batteries were built from compressed powders of 
Ko ^MoO^. In order for these materials to work as batteries 
the ions must pass from grain-boundary to grain-boundary. 
This implies that we do not even need single crystals (though 
no doubt epitaxy layers would be ideal) for preliminary 
experiments. Though in this situation we would not expect 
perfect control of the patterns of ions. Yet even with this we 
can exploit this phenomenon for computational purposes. The 
system, would act like a feedforward neural network [9 1 ] that 
can be programmed with a genetic algorithm. The program 
would be adjusting, essentially, tiny resistors. Figure IZSlshows 
a schematic diagram of our system and a prototype to test the 
basic ideas. 

Of course, if the material is polycrystal this will add un- 
controlled defects to our programmable surface. These crystal 
boundaries will result in faults in the computational network. 
However, as long as we take these faults into account with 
our design algorithms (a genetic algorithm) it should have 
little or no impact. Figure |25] shows a prototype system we 
built to test the basic idea of programmable resistor arrays. A 
film of K0.3M0O3 was made with 5 wt% polyethylene oxide 
as a binder and deposited on a glass slide with previously 
deposited gold electrodes with wires attached. The wires were 
connected to a simple screw terminal for strain relief. The 
two end terminals were connected together and then connected 
to a dc field (20 volts). Terminals from opposite sides were 
then connected to a 4 192 A impedance analyzer and the ac 
impedance was measured as a function of frequency. Figure 
l26l shows these, very preUminary, results. 
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Fig. 24. Logistic- and circle-map adding machines from linear arrays (see Sinaha and Ditto, 1998) 




Fig. 25. Simple programmable resistor network from charge density wave material 
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Fig. 26. Results from prototype of the CDW device 



The results indicate that there was more than an order of 
magnitude change in the resistance with the dc field on versus 
the dc field off. No doubt if we used epitaxy films ( [92]; [93]) 
we could have better control of the ions. 



VI. Acoustic Modulation of Media for Optical 
Computing 

It is possible to use acoustic pulses, in particular standing 
waves, for modulating the density of soUds, Uquids, and gases. 
In this section we review work by Higginson et al [94], [95], 
that describe the use of acoustic standing waves to modulate 
the density of media. 

The interaction between sound and optics was first described 
in [96], and later elaborated on in [97]. Normally acousto- 
optic effects arise form traveling acoustic waves producing a 
periodic modulation in the refractive index of the medium. 
This effect is well known for surface acoustic wave devices 
and produces a grating. The traveling acoustic wave has a 
much lower velocity relative to the speed of sound so it 
effectively acts as a grating. The effect has been exploited for 
an acousto-optic television display [98]. This television appli- 
cation suggests other technologies such as optical switching 
and mask-less Uthography. These systems exploit a traveUng 
wave. Contrast this with the work in [94], [95] who have been 
making use of an acoustic standing wave. 
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transducer 



Fig. 27. Acusto-optic device for lensing, described by [94], [95] 

In order to exploit an acoustic standing wave for modulating 
light they started with cylindrical device containing fluid 
(see Figure |27] i. One starts the analysis by considering the 
time-independent component of a standing wave, which can 
be derived from the Tait equation for liquids [99]). After 
expanding by a second order Taylor series we get 

p-po = ^(P-Po) + ^(P-Po)^ + ... 

where P and p are the total pressure and density,Poj Po, co 
are the ambient pressure, density, and sound speed. 7 is an 
empirical constant. 

Since the pressure is adiabatic, and the flow is irrotational, in 
the cylindrical cavity, expanding the total pressure in a Taylor 
series and grouping the terms can estimate the pressure terms. 
Time-averaging eliminates the first-order terms and gives the 
time-invariant component of the density in the liquid. 

where p and u are first-order acoustic pressure and velocity, 
obtained by solving the linearized wave equation, and the 
angular brackets indicate time average. Since the geometry 
is circular the acoustic standing wave can be described by 

p{r, t) = AJo{kr) sin{at) 

So we get p2{r) = ^ [(2 - 7) Jo(fcr) - Jfikr)] 

where A is the amplitude of the sound pressure, a is the 
angular frequency, k is the dispersion-free wave number u) / cq, 
Ji are Bessel functions of the first kind, and r is the radial 
coordinate. The density is related to the refractive index n, 
by Lorenz-Lorenz equation. We can write the time averaged 
version of this as 



(n) = no + a [(2 - 7) Jo(fcr) - Jl{kr)\ 

A^Q(nl + 2f 

(y. 

24p - pqcIuq 

where Q is the molar refractivity of the material, and n is 
the refractive index at po- The calculated time-averaged index 
profile of a glycerin-filled circular cavity is shown in Figure 




n.On D.25 0.511 



rjcrm) 

Fig. 28. Calculated refractive index profile in glycerin for A=200 MPa and 
n 27r = 700kHz . The inset is the refractive index disturbance as a function 
of acoustic pressure. 



l28l Here we have used the values 

po = 1260% 

Co = 1904m s 

no = 1.4746 

Q = 2.23X10"W kg 

7 = 10.0 

n2TT^ 700kHz 

A = IbOMPa 

Rather than obtaining as smooth gradient in the refractive, 
one finds that the maximum is at the center and there a 
decreasing periodic modulation in the refractive index, as 
suggested in Figure|28] These refractive index modulations can 
be observed by deflection of a laser beam. Figure |29] shows 
the results for this experiment. 

Figure |28] shows the calculated refractive index profile for 
glycerin with an assumed pressure of 200 MPa in the center. 
The refractive index change is on the order of 0.1%, but 
as shown in Figure |29] these seemingly small changes in 
the refractive index are large enough to deflect the beam by 
over a millimeter Based on these deflections, the pressure is 
determined to be about 50 MPa. Further, as shown in both 
Figure |28] and |29] there are radial modulations in the refractive 
index and these of course result in modulations of the light as 
shown in the beam-profile of Figure |29l 

These circular refractive index modulations in the ring 
produce an overall global behavior known as an axicon lens. 
Axicon lenses are used to produce optical beams with a 
Bessel and Mathieu function profile (Figure |29] l. An axicon 
fit is shown in Figure |28] These beams are said to be self- 
healing and diffraction-free. The intensity of the profile does 
not diverge diffractively as a function of distance from the 
lens, and diffraction around objects is eliminated. The beam 
is self-healing because the intensity pattern re-establishes itself 
some distance away from an obstruction. 
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Fig. 29. Predicted and actual beam deflection as a function of position. A 
= 50Mpa, at 412kHz. The photo on the right shows the overall beam profile 
at 30cm from the lens. This is a Bessel beam profile. 




Fig. 30. The figure on the left shows an acoustic driven lens with sectioned 
transducers that can be diiven independently. The top two photos show two 
of the many Hght patterns created and the lower two photos show two of the 
many standing wave pressure modulations possible in the fluid lens. 



More recent work has used multitransducers and sectioned 
ring transducers to provide more degrees of freedom in the 
fluid modulations. Figure [30] shows a photo of the device and 
several other photos. The bottom two photos were made by 
injecting 6 micron polystyrene spheres into the cell filled with 
water The individual transducers were powered up separately, 
and by a process know as acoustophoresis the micro-scale 
particles will collect at the pressure nodes in the acoustic 
standing wave. A mm-scale has been placed on the lens for 
reference purposes. 

The upper two photos are examples of the light patterns 
generated on the output side of the lens while it is powered 
up at different frequencies. It is likely that an infinite number 
of possible image patterns can be generated from the equally 
infinite number of particle patterns in the fluid. It is possible, 
to drive each of the transducers at a different frequency, to 
amplitude or frequency modulate the drive frequency at any 
one of the transducers, and to drive the transducers in some 
predefined sequence to produce desired optical images at the 
output. Naturally, with this many degrees of freedom, we 
expect to use a genetic algorithm to evolve the transducer drive 
frequencies and patterns. 

VII. Conclusions 

The results presented in this paper demonstrate that we 
can use bulk matter for computations almost as if it was 
programmable matter With a programmable "sea of gates" in 
silicon we evolved an analogue oscillating circuit. In a liquid 
crystal substrate we evolved a tone discriminator circuit. In a 



magnetic dot substrate we evolved (in numerical simulations) a 
pattern matching system. In more speculative work, we have 
modelled a programmable Fermi surface where we manip- 
ulated the charge density wave in a crystal lattice. Finally, 
we outlined some laboratory work in using acoustic pulses to 
modulate bulk matter to affect the refractive index of light for 
using the "constructed" optical device in optical computing. 

In each of the examples the computation comes about 
through local and nearest neighbour interactions, and the data 
and the program for the computations are fed directly into 
the bulk matter through the edges. If we had detailed models 
of the nano-scale physics, the computational burden would be 
significantly beyond foreseeable computer hardware. This can 
be circumvented through evolutionary programming. There are 
likely to be many more computational substrates that can be 
exploited through evolutionary engineering. 
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